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(§S) Signaltrennungsverfahren und -einrichtung fur nichtlineare Mischungen unbekannter Signale 

(g) Die Erfindung ist ein Verfahren und eine Einrichtung zur 
Separierung ernes Signalgemisches, bestehend aus der 
nichtlinearen Oberlagerung von M unbekannten Signalen, 
dadurch gekennzeichnet, das N (N & M) unterschiedliche 
Mischungen der M Quellsignale inklusive eines eventuellen 
Storsignals einer Signalauswerteeinrichtung zugefuhrt wer- 
den, die durch eine statistische Analyse der Signale die 
nichtlinearen Ubertragungsfaktoren bestimmt und mit die- 
sen errechneten Faktoren die Mischung ruckgangig macht, 
so da& die N Ausgange der Signaltrennungseinrichtung die 
M Quellsignale ohne Oberlagerungen enthalten. Dadurch 
wird erstmals die Behandlung nichtlinearer Gemische mog- 
Itch, wobei nichtlinear bedeutet, daS die Quellsignale durch 
ein unbekanntes nichtlineares System gemischt werden. Das 
unbekannte System wird durch eine Volterra-Reihe be- 
■ schrieben, und die Signaltrennungseinrichtung bestimmt die 

>Koeffizienten der Volterra-Reihe. Mit derer Kenntnis ist eine 
Entmischung moglich, auSerdem konnen die Koeffizienten in 
% weiteren Analysen zur Orts- oder Geschwindigkeitsbestim- 
mung der Signalquellen benutzt werden. 
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Beschreibung 



1. Technisches Gebiet 

In der MeQtechnik besteht fast immer das Problem, daQ bei einer Messung nicht nur Signale einer einzelnen 
Quelle, sondern oft nur Mischungen aus den Signalen mehrerer Quellen durch die Detektoren aufgenommen 
werden konnen. 

2. Stand derTechnik 

Eine Signaltrennung mit mdglichst wenig a-priori- Information uber die Signale wurde von Herault und Jutten 
in den folgenden Veroffentlichungen behandelt: C. Jutten and J. Herault, "Blind separation of sources, part I: An 
adaptive algorithm based on neuromatic architecture", Signal Processing, vol 24, No. 1 pp 1 — 10, July 1991; 
R Common, C Jutten, J. Herault, "Blind Separation of sources, part II: Problems statement", Signal Processing, 
vol 24, No. 1, pp 11—20, July 1991, E. Sorouchyari. C. Jutten, J. Herault, "Blind Separation of Sources, part III: 
Stability analysis, ", Signal Processing, vol 24, No. 1, pp 21-29, July 1991. Dabei wird das Problem von mit 
konstanten Faktoren linear iiberlagerten, zeitlich unverzogerten, unabhangigen Signalen beschrieben. Wie sich 
herausstellt, funktioniert das Verfahren nicht mit GauBverteilten Signalen, das sind zum Beispiel Sprache oder 
Musik. 

In der Veroffentlichung von Molgedey und Schuster (Molgedey L und Schuster H. G. 1994: Separation of a 
Mixture of Independent Signals using Time Delayed Corelations, Phys. Rev. Lett 72, 3634—3638.) wird die 
Trennung fur unkorrelierte Signale, die zeitlich unverzogert linear uberlagert sind, behandelt, die auch fur 
GauB'sche Quellen benutzt werden kann. AuBerdera wird ein Verfahren zur Trennung von Mischung skizziert, 
bei denen die Quellsignale (zeitlich unverzogert) multipliziert sind 

Gilles Burell (Burell, Gilles 1992: Blind Separation of Sources: A Nonlinear Neural Algorithm, Neural Net- 
works, Vol.5, pp 937—.) schlagt einen weiteren Ansatz fur die Trennung von Gemischen aus unabhangigen 
Signalen vor, in dem ein neuronales Netz speziell fur eine bekannte Art der nichtlinearen, zeitlich unverzogerten 
Mischungen konstruiert wird. Das Verfahren zum Finden der richtigen Parameter fur die Entmischung ist dem 
Backpropagationalgorithrnus ahnlich. 

In dem European Patent application number 93630032.6, date of filing 10.04. 1992 by the inventors Ehud 
Weinstein and Daniel Yellin wird eine Signaltrennung fur den Fall einer durch lineare Filter entstandenen 
Mischung unabhangiger und unkorrelierter Quellsignale beschrieben, und es werden zwei Berechnungsmdglich- 
keiten fur die gesuchten Filter gegeben: 1. Eine iterative Losung zweier "standard least square" Algorithmen und 
2. ein rekursiver Algorithmus, der im Prinzip dem ersten aquivalent ist, aber weniger Rechenkapazitat benotigt, 
dafur aber auch nicht exakt ist Damit wird der von Piatt und Faggin in (J. C Piatt and F. Faggin: Networks for 
the Separation of Sources that are Superimposed and Delayed, ed. John E. Moody "Advances in Neural Informa- 
tion Processing" Vol. 4, Morgan- Kaufman, San Mateo (1992).) vorgeschlagenen Algorithmus begrundet und 
erweitert. In keiner uns bekannten Veroffentlichung wird das Problem einer a priori unbekannten nichtlinearen 
Mischung (Mischung durch ein beliebiges nichtlineares System) angeschnitten oder gar geldst 

3. Zugrundeliegendes Problem 

Die Erfindung betrifft ein Verfahren sowie eine Einrichtung zur Signaltrennung eines Signalgemisches der 
Form 

i=i J0 

M rco /-co 

+ / / i72.Jt/( w i,"2)a:jL-(f -u x )x/(£ - u 2 )du l du 2 

+ ... 

M roo poo 
+ ZZ J '"J 9n % m x m„(«l,---i U n) 

das aus M unbekannten Quellsignalen x t (t) . . x M (t) besteht, die mit unbekannten Gewichtsfaktoren g 5 ^ ji.-, in ^ 
der angegebenen Weise gemischt sind. Ein Schema dazu stellt Fig. Nr. 1 dar. 1 steht fQr den unbekannten Teil, 
bestehend aus den Quellprozessen, deren Signale durch ein nichtlineares System G 2 gemischt werden. 3 zeigt 
den bekannten Bereich an. In 4 werden aus einer Analyse der MeBdaten die Quellsignale entmischt und stehen 
ais Ausgange u ( und u 2 zur Verfiigung. Das Verfahren kommt fast ohne Annahme uber die Quellsignale und das 
nichtiineare Mischungssystem aus, bis auf die Bedingung, daB eine statistische Analyse der empfangenen Mi- 
schungen aus den Quellsignalen ys (t) genugend Information fur eine Trennung liefert. Beispielsweise wird die 
zeitlich unverzogerte Oberlagerung zweier GauBverteilter weiBer Rauschsignale ni(t) und n^t) mit unterscnied- 
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lichen konstanten Faktoren C12 und c 2 i in der Form 
yi(t) = ni(t)+c 12 n2(t) 

y 2 (t) = n 2 {t) + c,2n l (t) (1) s 

nicht mogiich sein, da (yt(t)y2{t + T)) = 0 fur x*0 und alle Korrelationen der Form (yMOy^t + t)) = 0 fur 
m + n>2. . _ 

Die Beschrankungen werden im Verlauf dieser Beschreibung klarer werden, aber da eine Veranderung der 
g'n, mi mn durch eine geanderte Auf nahmetechnik radglich ist und da man es in den meisten Fallen rait der 
Trennung von Signalen, die Information enthalten, d. h. bei denen eine statistische Analyse hoherer Ordnung io 
Ergebnisse liefert, wie z. B. eine Trennung von Mischungen aus Sprache und Musik, zu tun hat, werden die 
Beschrankungen meistens ausraumbar bzw. vernachlassigbar sein. 

4. Erfindnung 

15 

Eine beliebige nichdineare Mischung kann durch die folgende Form beschrieben werden: 
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wobei die M Quellsignale *i(t) . . . xm(0 unbekannt und mit unbekannten Gewichtsfaktoren g 5 ^ j t# .. ,j n in der 
angegebenen Weise gemischt sind und wobei teR ist Dies ist die Form einer Vol terra- Reihe. Entsprechend 
einem Theorem von Frechet kann jedes stetige Funktional mit einer beliebigen Genauigkeit durch eine Volter- 35 
ra- Reihe approximiert werden. Zu dieser Reihe gibt es auch ein zeitdiskontinuierliches Analogon, das eine 
Anwendung der Erfindung auf digitale Systeme zulaBt 

Im folgenden wird nur der Fall betrachtet, in dem xj(t), ys(t), teR gilt Eine Formulierung mit hoherdimensio- 
nalen Variablen ist mogiich, aber der eindimensionale Fall enthalt alle wesentlichen Verf ahrensschritte. 

Das in dieser Erfindung beschriebene Verfahren besteht im wesentlichen aus zwei Schritten: Erstens werden 40 
die durch den wahlbaren Grad der Nichtlinearitat bei der Mischung eindeutig bestimmten nichtlinearen Glei- 
chungen fQr ein gleitendes Zeitfenster gelost und die Losungen werden uber die Zeit gemittelt, und zweitens 
wird ein aus einer geniigend groBen Anzahl von unterschiedlichen Kumulanten der geschatzten Ausgangssigna- 
le uj(t) gebildetes Potential minimiert, wobei die zur Berechnung des Potentials ndtigen Werte aus einem 
gleitenden Zeitfenster einer wahlbaren Lange stammerL 

Zu erstens: 

Eine mehrdimensionale Fouriertransformation von 2 ist mogiich, und man erhalt: 



M 50 

+ £ G'n.m, («^i,..-,«--)A' mi (a; 1 )...A' mil {u: B ). 

mi »in = l 55 

Genauso ist eine zeitvariable Fouriertransformation mdglich, in der jeweils ein zeitlich begrenzter Ausschnitt 
der GroBe TeR aus dem Gesamtverlauf der Signale herausgeschnitten und fouriertransformiert wird. Dann 
erhalt man: 

60 

yaz.t) = i:<r^(u, l )A' i (^,r) + ... 

+ £ <3L« mB (^,...,w„)A' mi (w,,r)...A' mn (w n ,r). 
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Dabei ist vorausgesetzt, daB sich das System G so langsam andert, daB diese Anderung bei der Berechnung der 
gesuchten Mischungsfaktoren vernachlassigt werden kann. 

Urn die Bedingungen zu formulieren, die die Quellsignale erfullen mussen, damit das Verfahren im ersten 
Schritt angewendet werden kann, werden nun folgende Variablen definiert: 

5 ij,k.No,N,^J 2 eN 
coij€R 

(i,,..,iNo)€{peN|l^p<M} N o 
to Oi— jM,)6N N i 
(k,...,k N2 )€N N 2 
(oSjjp.^OiNojNOeRNoNl 

( *iuiM. • • - TiNojNik N2 ) e RNoNiN 2 

15 Die Quellensignale mussen die folgenden Bedingungen erfullen fur einen genugend groBen Bereich der oben 
definierten Variablen: 



20 



cwm( (3) 

AVj (^'m jni » ^ ~* r »i y.v,^i )*•••- A (wfjj^ , T — TV, Jjif| ) , 

A"iV 0 (^l^pii » T r T U^JtMi )i ■ • - 1 A% 0 (o><„ oil ,T - Ti NoJl ^ ) t 

A'/zw^C^fVoAv, ' T ~~ r »v 0 J.v,,*ri )••*• • Ai^^f^/^t* r ^v 6 J.v, ) 
) = 0. 

35 auBer falls 

i, = i 2 . . . = i No 
und 
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C0i!j| = ... » OiNChjNl 



gleichzeitig gilt 

Nun geht es darum, die in der obigen Gleichung stehenden unbekannten GroBen zu berechnen. Man bendtigt 
45 dieselbe Anzahl an Gleichungen wie an unbekannte GroBen. Daher soil mit obiger Gleichung fur ausreichend 
unterschiedlich viele Verzogerungszeiten t 0 bis te(E€N) ein Gleichungssystem aufgestellt werden, d. h,: 



50 



(Ys 0 {Z,T)Y Si (<Z^T + T 0 )) = ... 

{Y^T)Y^T + r E )) = ... (4) 

55 oder entsprechende Gleichungen, die sich durch das Bilden der Erwartungswerte von beliebigen Produkten der 
zeitvariablen Fouriertransformierten der empfangenen Signale Ys(a>i, . . ., G> n , T) und der zeitverzogerten zeitva- 
riablen Fouriertransformierten der empfangenen Signale Y s (g>i, . . , o>n, T— Xj) ergebea Dadurch erhalt man ein 
nichtlineares Gleichungssystem, das numerisch gelost werden kann. Naturlich kann man durch die Berucksichti- 
gung weiterer Werte fur die Verzogerungszeiten xj mehr Gleichungen als unbedingt notig aufstellen und dann 
60 die uberbestimmten Losungen einer Mittelwertbildung unterziehen. 

An einem Beispiel soli nun die Verfahrensweise dieses ersten Schrittes gezeigt werden. 
Das betrachtete System sei derart nichtlinear, daB gilt: 



65 
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Y^.u^T) = Xi(«i.r) + 

G*!. 2 (^)X,(u.-,.r) 

Gl, .,(»,•, , (c,-» - T)X,(^, T) 
Cl^Wu^XiWu T)X X {^, T). 

Dann ergeben sich durch das Bilden von Kreuzkorrelationen die folgenden Gleichungen (0<j<E) fOr 
ct>2 und o>i •Aon: 

{Y|(oi,<i)2,T)Y2(aJi.a)2,T+Tj)> = 
G 2 i.i(coiXXi(coi,T)Xi((Oi,T+T j )>+ 
G> i^(o)i YX2((i)i, T)X2(t0i, T + Tj)> + 

G' 2 /^o>V,co2)G^ 1 ^a),,o> 2 XX 1 (co,,T)X 1 ((0,.T+T j )XXKo>2,W 
G«2^(©i,co2)G 2 2Ai(o>,,<02XX 

(Yi(cai,o»2,T)Yi(fi)i,tD2,T+Tj)) = 

<X,(oa,,T)X,(o)i,T+T)> + 

(G't^)) 2 (X2{to,,T)X2(o>i,T+T j )>+ 

(G't,!^©,, asHMO" T)Xi(<oi,T + Tj)XX2<<02, T + Tj)) + 

{G'2Al(fi>l. ©2)) 2 (X2(<01, T)X2(ft)l,T + T j )XX,(O2.T)X,(0>2, T+ Tj)); 

(X2(<a,,T)X2{fi)i,T+Tj)>+ 
(G J u(©,))2(X,(o>i.T)X,((D,,T+Tj)) + 

(G 2 !^©,, ©zJ^.Coa,, T)X,(©,,T + T j )XX 2 (o)2,T)X 2 (o>2.T + tj)>+ 
(G 2 2Ai(©i.<^)) 2 ^2(fi>t.T)X2(cai,T+Tj)XX,((02.'T)X,(a>2,T+Tj)>. 

Entsprechend ergeben sich 3 weitere Gleichungen, wenn on und <o 2 vertauscht werden. 

(Y,(<a 2 ,<0i.T)Y2(<ai,c82.T + Tj)> = 
G 2 t,i(o)2XX I (o)2.T)X,(fi> 2 . T +Tj))+ 

G'i;K©2XX2(<B2,T)X2(G>2,T + Tj))+ _ , ~ « 

G'2.u(<02, to.^i^, a),XX,(6>2,T)X 1 (©2,T+ ^XXKffiuTJXaCaM.T+Tj)) + 
G'2Ai(ffl2, aiJG^itffl* ©iXX2(o> 2 , T)X2(a)2, T + Tj)XX,(o>,, T)X,(cd,. T + Tj)>; 

(Yi(<fl2,a)i,T)yi(oi>2,<o,,T+Tj)) = 

{X,(d)2.T)X 1 (0>2,T + Tj)> + 
(G'l^2)) 2 (X2( 0 >2,T)X2(0>2,T + Tj)> + 

(G'zi^oz, <Di)ftXi(«»» T)X,((02.T+ Tj)XX2(a>,,T)X2(G>,, T+ *j)>+ 
(G'2Ai(w2, ».))HX«(a>2. T)X 2 (<02,T+T j )XX,(©wT)X 1 (o>,.T+Tj)); 

(Y2{0)2. Q>|, T)Y2(0>2, ©1, T + Tj)) = 
(X2(0)2.T)X2((a2.T + Tj)) + 

(G 2 l,l(0) 2 )) 2 {X,(<D2,T)X 1 (<02,T + Tj))+ 

(G*2.1,2(C02. ©.WW**, T)X,(<02, T + TjJXX^OJ,. W©,. T + Tj)) + 
(G 2 2Al( m 2. O),)) 2 ^^, T)X2(<02. T + Tj)XX,(0>,, T)X,((D,, T + Tj)). 

Fur ein Paar (to,, ©2) mit ©, # <o 2 ergeben sich 6E Gleichungen fur 12 + 4E unbekannte GroBen. 
InTFall, daB L^eR mit o>, =0* ist. ergeben sich durch das Bilden von Kreuzkorrelationen die folgenden 
Gleichungen (OS j S E): 

(Yi(ffli.o>i,T)Y2(oi>i,<D,,T + Tj)) = 
(G 2 i.i(ft>,XXi(o)i,T)X 1 («i)t,T+Tj)) + 
(G I i^(a)i)(X2(a)i,T)X2(oi)i,T+Tj)) + 
fe« 4 1k»Y.«oi)G^i*»i,a>^ 

G' 2 Ai(©i.©i)G 2 2Ai(tOi.<OiXX2(to,.T)X2((0,.T+Tj)XX,(© 1 ,T)X,(co,,T + ^^^ 

G'2 >1 ^<O,.<O,)G 2 2il(ffi t> tD2XX2(©.,T)X2(a> 1 ,T + Tj)^ 1 (0 i| ,^X,t0 

G^K^i. ^i)° 2 ii^ ffl «' w 'XX2(o>i, T)X 2 (&) 1( T+ TjJXX.Cco,, T)X,(to,.T + Tj)); 
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(Y,(£Oi,ai.T)Y,(G>i,o>i.T+Tj)) = 

(X,(tDi,T)X,(a>,,T+Tj)> + 

(G 1 i^(o,)) 2 (X 2 (co,.T)X 2 ((D,.T+Tj)>+ 

(G'2.ii(o>i.co,)) 2 (X,(a>,,T)X 1 (o),,T+T j ))(X2(o)i,T)X 2 (co,.T+Tj))+ 
(G , 2Ai(wi,o),)) 2 pC2((o,.T)X 2 (o),,T+Tj)XX,(o)i,T)X,(ti>i,T+Tj))+ 
G , 2.j^©i,o) l )G , 2Ai(© l ,w,XX,(co,,T)X,(o>,.T + T j )XX 2 (a),,T)X2(o> l ,T+Tj)) + 
G'2.i^(o),, co,)G l 2Ai(o)i , o,XX2(a), . TJX^cowT + Tj)XX,(a>,, T)X,(<0,. T + Tj)>; 

{Y 2 (a>i,o>i,T)Y^o>,,(Oi,T+Tj)) = 

{X 2 ((B|.T)X2( tI >,,T + T j )) + 

(^i.i^OXXiCoii. T)X,(a>,,T + tj)> + 

(G 2 2.ii(co, > o>,)) 2 (X,((o,,T)X l (o),,T+Tj)XX2(o t ,T)X2(o)i,T+Tj)>+ 
(G 2 2Ai(«>i. o>i)) J <X2(co,. TJXKoi. T + T^XX.Cto, , T)X,(a>,, T + tj)) + 
G^i^ca,, <o,)G 2 2Ai(©i . <o,XX,(o),, T)X,(oi,T + Tj)XX 2 (a>,, T)X 2 (<o,, T + Tj)> + 
G 2 2,i^Oi,o)i)G 2 2Ai(<»i. m iXXKwi»^X 2 ((ii|,T + Tj)XX|(a)i,T)X|(ti>i,T+Tj)). 

Fur ein Paar (o>i, o> 2 ) mit a>, = cd 2 ergeben sich 3E Gleichungen f ur 6 + 2E unbekannte GroBen. 
Beim Aufstellen dieses Gleichungssystems wurde davon Gebrauch gemacht, daB die Quellsignale die folgen- 
den Bedingungen entsprechend Gleichung (3) erf Alien: 



cum(Xi (uj, T), Xi{u>i. T — 


r)) 




0 


cum(.Yi(u.* 1> r),.Vi(w l ,r — 


t),X 2 (us 7 ,T — 


r)) 
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cum(.V 2 (w, , T), A'j (w, , T — 


t), A' 2 (w 2t r — 
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cum(A'i {t^,T), A" 2 (u;i , T - 


r), Ai(u> 2 , r - 


r» 
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cum(A' 2 (fa>i , 7*), A' 2 (u;i , T — 




r)) 




0 


cum(A'[(wi,r), A' 2 (ua, T), A^u;, , T - 


r), A' 2 (a^,T - 


r)) 




0 


cum(.Yi(w, , T), A' 2 (u; 2 , T), A' 2 (w, , T - 




r)) 




0 


cum( A' 2 ("i , T), A'i (a.v,, T), A' 2 (a,, , T - 


t),Xi[u} 2 -,T — 


t)). 




0 


ctx»n(A' 2 (u> 1 ,r),A'i(u.'2, D, -Y,(u>i, T - 


r),A'2(wj,r — 


r)) 




0 



fur alle a>i, 02, t e R, und 

cum(Xi(coi,T),Xi(a> 2 ,T-T)) - 0 
cum(X 2 (a>i ( T) > X2((o 2 ,T-T)) = 0 

fur alle coi o>2, x e R mit ©i gfc 0)2. 
Da mit gilt dann: 

pC^couTJX^TJX^WtT-tJX^T-T)) = 
(Xi(<Oi,T)Xi(a>i,T— tJXX^ T)X2(a)2, T- t)>; 
(X^OuTJX^o^TJX^couT-tJX^^T-t)) « 
(Xj(©i, TpC^ffit. T-t)XX,((o 2 , T)Xi(<02, T-t)>; 
(X^cD^TOX^caL-DX.^i.T-TjX^tOwT-x)) = 
(X2(g>i, T)X2<Ci>i ( T- T)XX,(co lf T)Xi(oi, T- 1)>; 

fur alle ti>i, o>2 € R. AuBerdem gilt: 

(X^cOuTJXKcoj.'DX^^T-TiXKwwT-T)) - 0; 
(Xa(©i f T)Xi(©2. T)X|(a>i,T-t)X2(o)i,T-T)) = 0; 

fiir alle a>t, (02 e R mit cot ^ a>2. 

Aufgrund von fehlerhaften Aufnahmen (zusatzliches Rauschen) bzw. Fensterungseffekten durch die zeitvan- 
ablen Fouriertransformation muB diese Berechnung fur ein in geeignet groBen Schritten uber die Sign ale 
bewegtes Fenster wiederhoh werden. Die entsprechenden Ergebnisse werden dann gemittelt, und die gemittel- 
ten Werte nahern sich so den wahren Werten an. Diese Annaherung kann aber viele Daten benotigen. Eine 
groBe Menge an Daten kann zu Problemen bei einer Realisierung fuhren, daher wird diesem ersten Schntt ein 
zweiter nachgestellt 

Zuzweitens: . 
Dieser zweite Schritt ist eine Anlehnung an die Methoden, die in der Signalverarbeitung mit neuronalen Netzen 
angewendet werden. Man konstruiert sich eine Kostenfunktion, die minimiert wird. Wenn das globale Minimum 
erreicht ist, hat man die optimalen Werte, in diesem Fall der Obertragungsfaktoren, gefunden. Fur die Mmimie- 



DE 195 31 388 CI 



rung einer Kostenfunktion stehen beim heutigen Stand der Technik eine groBe Anzahl von Methoden zur 
Verfugung, z. B. Gradientenabstieg oder Monte-Carlo-Simulationen. Ein Nachteil aller Verfahren ist aber, daB 
sie in lokalen Minima der Kostenfunktion steckenbleiben konnen; der Vorteil unserer Erfindung ist, daB durch 
die Startwerte aus dem ersten Schritt die Minimierung schon in der Nahe des globaien Minimums beginnen 
kann, denn bei dem Ende der Berechnungen des ersten Schrittes liegen fur die Obertragungsfaktoren geschatzte 5 
Werte vor, die sich nur wenig von den gesuchten Werten fur die wahren Obertragungsfaktoren unterscheiden 
und die fur eine bessere Trennung weiter optimiert werden. 

Mit der Methode der Potentialrainimierung werden Effekte von zusatzlichen Rauchprozessen verringert. 

Urn die Bedingungen zu formulieren, die die Quellsignale erfuJlen mussen, damit das Verfahren im zweiten 
Schritt angewendet werden kann, werden nun folgende Variablen definiert: to 

i,j,No,Nt€N (5) 
(i I( ...i No )e{p€N|l<p<M} N o 

ai,...jN,)eN N t is 
(r llil( ..,r lNO jNi)eR N o N i 



Die Quellensignale raussen die folgenden Bedingungen erfuJlen fur einen genugend groBen Bereich der oben 
definierten Variablen: 

cum( (6) 

X M U ~~ r i*i Ji )i • * - ? X h ~~ T M Jfti )> 

) = o ; 



auBer falls 



= >N 0 



20 



25 



30 



gilt. 35 

Fur den hier betrachteten Fall, die richtigen Obertragungsfaktoren zu linden, bilden wir eine Kostenfunktion 
derart, daB die Ausgange eines adaptiven Filters ui(t) dieselben Bedingungen erfiillen mussen wie die Quellsigna- 
le x;(t), d. h. ftir eine genugend groBe Anzahi von Werten der in 5 definierten Variablen werden die Quadrate der 
Kumulanten aus 6, dabei werden die Quellsignale x durch die Ausgange des Filters u ersetzt, aufsummiert, so daB 
eine Kostenfunktion entsteht, die minimiert wird Optimiert werden dabei die Werte der Obertragungsfaktoren. 40 
Nach jedem Optimierungsschritt mussen die Werte fur die Ui(t) neu berechnet werden. Diese Berechnung kann 
im Zeit- oder Frequenzbereich erfolgen. 

Zum Beispiel werden fur den Fall, in dem folgende Mischung vorliegt: 

Yi((i>) = X l (co) + Gi 2 (co)X2(co)+ F^coJX^oJX^o)); 45 

Y 2 (0)) = X2(Q)) + G2l(CD)Xi(<i))+K(C0)Xi(ci>)X2(C)) 

die Ausgange mit den geschatzten Werten fur die Obertragungsfaktoren G*i& G*2i und K* berechnet durch 
Aufldsen des Gleichungssystems: 

50 

Y,(<o) = Ui(6>) + G*i2(co)U2(a)) + K # (©)Ui((o)U2(a)); 
Y2(a>) = U^)^G*2i(co)Ui((o)^K\<^)Vi((o)U^) 

nach Ui(o>) und U^co). 

Die durch das adaptive Filter zu rainimierende Kostenfunktion V kann mit den Ausgangen Ui(t) und U2<t) 55 
gebildet werden: 



60 



[cum-iu^tlu^i + r),u 2 (t + r))] 2 + 
[ciim(u 2 (t),u 2 (t + r),ui(t + r))] 2 + 

[cu77*( Wl (0, « 2 (0> + T )* u 2(t + r))] 2 -h 65 
[cum(u 2 (0,ui(0i«2(* + r),u,(< + r))] 2 , 



7 



DE 195 31 388 CI 



wobei B eine endliche Menge reeller Zahlen ist. 

5. Gewerbliche Anwendbarkeit 

Die durchzufuhrenden Berechnungen bendtigen eine hohe Rechenkapazitat, die aber schon heute mit Mikro- 
prozessoren bereitsteht Eine in parallelen Schritten ausgefuhrte Verwendung des angegebene Verfahrens 
eroffnet auch fur haherfrequente Probleme ein Echtzeitrealisierung. Im ofnine Betneb laBt sich das Verfahren 
auf frei programmierbaren Mikroprozessoren implementieren. . 

Anwendungsmdglichkeiten liegen in einer Signaltrennung als Vorverarbeitungsschntt fur eine automatiscne 
Spracherkennung, bei der Konstruktion technischer Hdrhilfen, sowie in alien anderen akustischen Messungen, 
einschlieBlich die der Unterwasserakustik. . 

AuBerdem sind medizinische Anwendungsmoglichkeiten bei Mehrkanalmessungen wie EEG Oder MbO 
(Fig. Nr. 2) moglich, sowie bei Antennenfeldern fur die Aufnahme elektromagnetischer Strahlung. Wie t in der 
Fig. Nr. 2 dargestellt, werden magnetische Dipole im Gehirn 2 durch SQUIDs 3 gemessen, und die Signale 
werden einer Auswerteeinheit 1 zugefuhrt. . 

Diesen und alien weiteren denkbaren Anwendungsmoglichkeiten ist gemeinsam, daB sie der bisher linearen 
Theorie verschlossen bleiben, weil eine lineare Signaltrennung zu unbefriedigenden Ergebnissen fuhrt, denn 
aufgrund der Detektortechniken haben fast alle Empf anger eine nichtlineare Kennlinie. Weiterhin sind bei einer 
Formulierung mit mehrdimensionalen Variablen Anwendungen in der Optik bzw. Bitdverarbeitung moglich, in 
der auch bei der Aufnahme nichtlineare Effekte auftreten. ^ 

Daruber hinaus konnen die berechneten Werte auch zu einer Analyse des mischenden Systems benutzt 
werden, z. B. zur Quellenlokalisierung oder Bewegungsdetektion. 

6. Vorteilhafte Auswirkungen 

Neben den schon in dem Abschnitt Gewerbliche Anwendbarkeit genannten moglichen Erneuerungen gibt es 
noch eine weitere Verbesserung: Um mit der bisher linearen Methode zu arbeiten, werden Empfindhchkeitsbe- 
reiche der Detektoren auf ihren linearen Bereich verkleinert Dies ist bei der nichtlinearen Methode nicht notig. 

7. Ausfuhrung 

Stellvertretend fQr alle anderen Probleme wird der Fall derTrennung akustischer Signale behandelt Entspre- 
chend Fig. Nr. 1 stellt der mit 1 bezeichnete Block dann einen unbekannten ProzeB bestehend aus zwei Schall- 
quellen dar, deren Signale durch den Raum zu zwei Mikrofonen gelangen, deren aufgenommene Signale 
aufgezeichnet werden und die Grundlage der der Erfindung entsprechenden Berechnungen bilden. Der Block, 
der mit 3 gekennzeichnet ist, enthalt alle bekannten GroBen. Das System G ist nichtlinear, weil zura einen die 
Raumubertragung nichtlinear sein kann, zum anderen die Mikrofone mit ihren Verstarkern eine nichtlineare 
Kennlinie besitzen. Die berechneten Signale u t und u 2 sind den Quellsignalen Xi und x 2 ahnlich, so daB ein 
automatisches Spracherkennungssystem mit der aus der Mischung herausgetrennten Sprache funktiomert. Dies 
ist in der Fig. Nr. 3 sichtbar. Dort sind in den oberen beiden Kasten die Mischsignale von Sprache und Musik 
aufgezeigt, wie sie sich bei einer Mikrofonaufnahme in einem kJeinen Raum ergeben- Im untersten Kasten ist das 
aus dem Gemisch herausgetrennte Sprachsignal am Ausgang der Trennungseinrichtung, dessen Berechnung auf 
einem der Erfindung entsprechenden Verfahren basiert, dargestellt; mit dem untersten Signal ist eine fehlerlose 
automatische Erkennung mit einem kommerziellen automatischen Spracherkenner moglich, wahrend derselbe 
Erkenner bei den oberen beiden Signalen jeweils drei Fehler liefert 

8. Figurenbeschreibung 

Fig. Nr 1:1 steht fur den unbekannten Bereich, bestehend aus den Quellprozessen, deren Signale durch ein 
nichtlineares System G 2 gemischt werden. 3 zeigt den bekannten Bereich an. In 4 werden aus einer Analyse der 
MeBdaten die Quellsignale entmischt und stehen als Ausgange Ui und u 2 zur Verfugung. . 

Fig. Nr. 2: Magnetische Dipole im Gehirn 2 werden durch SQUIDs 3 gemessene, und die Signale werden einer 
Auswerteeinheit 1 zugefuhrt , _ r . . . . . 

Fie. Nr 3: In den oberen beiden Kasten sind die Mischsignale von Sprache und Musik aufgezeigt, wie sie sich 
bei einer Mikrofonaufnahme in einem kleinen Raum ergeben. Im untersten Kasten ist das aus dem Gemisch 
herausgetrennte Sprachsignal am Ausgang der Trennungseinrichtung, dessen Berechnung auf einem der Erfin. 
dung entsprechenden Verfahren basiert, dargestellt; mit dem untersten Signal ist eine fehlerlose automatische 
Erkennung mit einem kommerziellen automatischen Spracherkenner moglich, wahrend derselbe Erkenner bei 
den oberen beiden Signalen jeweils drei FehJer liefert 

Patentanspruche 
1. Verfahren zur Separierung eines Signalgemisches der Form 
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+ ... 

+ £ j~ - l~ 9 «™ m>l> ••••«») 10 

rn \ ...m n <=l 

Im, (* ~ " 1 ) • • • *m„ (* - «n)^l - - -du n , 

das aus M unbekannten Quellsignalen x^t)... xmO) besteht, die mit unbekannten Gewichtsfaktoren is 
g s n j!,..,j n in der angegebenen Weise gemischt sind, fur die Variable t soil gelten: teR, dadurch gekenn- 
zeichnet, daB ersten die durch den Grad der Nichtlinearitat bei der Mischung eindeutig bestimmten 
nichtlinearen Gleichungen, die durch das Bilden von einer genugend groBen Anzahl von unterschiedlichen 
Erwartungswerten von Produkten der zeitvariablen Fouriertransformierten der empfangenen Signale 
Y s (o) 1( . . , co m T) und der um xj zeitverzogerten zeitvariablen Fouriertransformierten der empfangenen 20 
Signale Y s (wi, . . <o n , T— xj) entstehen, fur ein gleitendes Zeitfenster geldst werden und die Losungen fiber 
die Zeit gemittelt werden und daB zweitens ein aus einer genugend groBen Anzahl von unterschiedlichen 
Kumulanten aus den geschatzten Ausgangssignalen Ui(t) und aus den um xj zeitverzogerten geschatzten 
Ausgangssignalen ui(t — xj) gebildetes Potential minimiert wird, wobei die zur Berechnung des Potentials 
notigen Werte aus einem gleitenden Zeitfenster einer wahlbaren GroBe stammen. 25 

2. Verf ahren nach Anspruch 1 , dadurch gekennzeichnet, daB fiber die Werte der Nichtlinearitat g irgendwel- 
che a-priori-Informationen bestehen, so daB die Berechnungen vereinfacht werden kdnnen. 

3. Verf ahren nach Anspruch 1, dadurch gekennzeichnet, daB M = 2 ist 

4. Verfahren nach Anspruch 1 bis 3, dadurch gekennzeichnet, daB fiber die Werte der Quellsignale x 
irgendwelche a-priori-Informationen bestehen, so daB die Berechnungen vereinfacht werden konnen. 30 

5. Einrichtungen, dadurch gekennzeichnet, daB sie die Verfahren der Anspriiche 1 bis 4 in digitaler oder 
analoger Weise fur eine Signal trennung benutzen. 

Hierzu 2 Seite(n) Zeichnungen 
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